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I. INTRODUCTION 

It is an old idea that the light scalar mesons a(980) 
and /(980) may be 4-quark bound states yj. The idea 
was more or less accepted in the mid-seventies but then it 
lost momentum, due to contradictory results that led the 
lowest-lying candidate members of a diquark-antidiquark 
nonet, a and k, to disappear from the Particle Data Ta- 
bles. 

As an alternative, the possibility was considered that 
a(980) and /(980) may be K — K bound states, kept to- 
gether by hadron exchange forces, the same that bind nu- 
clcons in the nuclei, color singlet remnants of the confin- 
ing color forces (hence the name K — K molecules Q used 
in this connection). If they are indeed K — K molecules, 
scalar mesons do not need to make a complete SU(3) mul- 
tiplet so that this idea would be consistent with the lack 
of evidence of light a and k. On the contrary, since the 
latter particles would in any case lie considerably higher 
than the respective thresholds, it would be very hard to 
consider either of them as a i - 7ror7r — K molecule. 
We see that the existence or absence of the light scalars 
is crucial in assessing the nature of a (980) and /(980). 
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From this point of view, the recent observations of 
(j(480) and k(800) in D non-leptonic decays at Fermi- 
lab H and in the 7T7t spectrum in <j> — > 7r°7r°7 at Fras- 
cati [J, have considerably reinforced the case of a full 
nonet with inverted spectrum, as expected for [qq] [qq] 
states and fully antisymmetric diquark ([qq] : color = 
3, flavor = 3, spin = 0). The isolated 1 = state is the 
lightest and it likes to decay in 7T7t; the heaviest parti- 
cles have I = 1,0 and like to decay in states containing 
strange quark pairs. 

In the late seventies, diquark-antidiquark mesons have 
been considered also from a different point of view, the 
so-called baryonium [j| . Baryonium resonances are called 
for by the extension to baryon-antibaryon scattering of 
the Harari-Rosner duality |f|, which is well obeyed by 
meson-baryon amplitudes. In the latter case, the ex- 
change of qq mesons in the t-channel is dual to (indeed it 
implies the existence of) non-exotic baryon resonances in 
the s-channel. Similarly, qq exchange in the t-channel of 
baryon-antibaryon scattering would give rise to 
states in the s-channel. 

In the constituent quark model, emphasis for the bind- 
ing of the diquark is on spin-spin forces, which may 
lead to strong attraction in the completely antisymmet- 
ric state. The baryonium picture adds a further ele- 
ment, the internal string structure associated with the 
confining, spin- independent, color forces jjj. QCD vac- 
uum restricts to a one-dimensional string the color lines 
of force emerging from each quark Q. In baryons, the 
strings from the three quarks join in a point to form a 
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gauge invariant color singlet and give rise to a V-shaped 
topological structure, with a quark sitting at each ends 
of the Y. Decay of baryon resonances is produced by 
the breaking of one of these strings, to give a qq me- 
son and a lighter baryon (e.g. A — > tt + N). In the 
same picture, duality implies the string structure of the 
diquark-antidiquark states to be that of an H, with the 
two quarks sitting on one side and the two antiquarks 
on the other side of the H . Topologically, an iJ-shaped 
state can be seen to arise from the fusion of two V-shaped 
objects, the baryon-antibaryon pair. Conversely, OZI al- 
lowed decays originating from the breaking of a string 
correspond to decays into either a lighter qqqq state plus 
a qq meson or into baryon-antibaryon. For the lightest 
scalar mesons of each flavor these channels are forbidden 
by energy conservation and we expect basically narrow 
states. The decay into meson-meson pairs has to proceed 
via the tunneling of the iJ-shaped configuration into two 
particles with quarks and antiquarks joined by a single 
string: S — ► (qq) + (qq). This picture is shown in [i| to 
give a reasonable description of the decay amplitudes of 
the lightest scalar mesons. 

If the lightest scalar mesons are diquark-antidiquark 
composites, it is natural to consider analogous states with 
one or more heavy constituents @ ^3] (see for an 
early proposal). 

The aim of the present paper is a study of diquark- 
antidiquark states with hidden or open charm of the 
form: [c^l[cq'] and [cgjfstf], q,q' = u,d. With respect 
to Ref. [9j, we add two new elements, the near spin- 
independence of heavy quark forces and isospin breaking 
from quark masses. We find some unexpected results and 
predictions, summarized in the following. 

For [cq] [cq'] states, the approximate spin-independence 
of heavy quark interactions |l2j . which is exact in the 
limit of infinite charm mass, implies spin one diquarks 
to form bound states if spin zero diquarks do so ("bad" 
and "good" diquarks, in Jaffe's terminology [13 )■ A rich 
spectrum is implied, with states with J — 0,1,2 and 
both natural and unnatural J PC . We describe the mass 
spectrum in terms of (i) the constituent diquark mass 
and (ii) spin-spin interactions. We derive the strength of 
the latter interactions from the known meson and baryon 
spectrum, where possible, or from educated guesses from 
one-gluon exchange, otherwise. 

We identify the V(3872) with the J PC = 
1 ++ state with the symmetric spin distribution 
[cq] s =i[cq}s=Q + [cq]s=o[cq]s=i (the charmonium assign- 
ment and its difficulties are described in Ref. Our 
assignment is consistent with the observed decays into 
charmonium plus vector mesons. It also implies a pure 
5 = 1 configuration for the c — c pair, thus complying 
with the selection rules derived in Ref. • We have one 
J PC = 2 ++ state at 3952 MeV that, within the accuracy 
of the model, could be identified with the AT (3940) seen 
in Belle data The scheme features two, J PC = ++ , 
states not yet identified. One, at 3830 MeV, could decay 
into D — D while the other is below the D — D thresh- 



old and should decay into rj c + ps mesons or multihadron 
states. Finally, there are two J PC = l" 1 levels, pre- 
dicted around 3760 MeV and 3880 MeV, also not yet 
seen. 

It is unclear to us if "bad" diquarks with light flavors 
can bind to [cq], let alone to a completely light-flavored 
antidiquark, or if the stronger repulsion in the S = 1 
state will suppress bound state formation. In the first 
case, an even richer spectrum is implied, due to the flavor 
symmetry of the light diquark, with many exotic states. 

We extend the previous calculation to the spectrum 
of [cq] [sq'] states, which can be computed on the basis 
of the same parameters. The resulting spectrum can ac- 
commodate the A (2632) claimed by the SELEX Collab- 
oration as well as two other states previously dis- 
covered, namely D s (2317) (with J p = 0+) and D s (2457) 
(J p = 1 + ) [l^, which could be at odds with a cs assigna- 
tion [23| • The latter hypothesis has been discussed in var- 
ious papers [2l]]. Also a molecular composition of these 
states has been taken in consideration, see for example 
Ref. ■ For a review and a more extended collection of 
references on this topic see e.g. [23). 

At the large momentum scales implied by the heavy 
quark, the strength of self-energy annihilation diagrams 
decreases. As a consequence, particle masses should be 
approximately diagonal with quark masses, even for the 
up and down quarks [lOl |24j. Neglecting annihilation 
diagrams, the neutral mass eigenstates coincide with: 

X u = [cu][cu]; X d = [cd][cd) (1) 

Deviations from this ideal situation is described by a mix- 
ing angle between X u and Xd- Considering the higher (h) 
and lower (1) eigenvalues, we predict: 

M(X h ) - M(Xi) = 2(m d - m„)/cos(20) = 

= (7±2)/cos(20) MeV (2) 

in terms of the up and down quark mass difference [2^ ■ 

Isospin is broken in the mass eigenstates and, conse- 
quently, in their strong decays. In particular, we expect 
this to be the case for Xh and Xi, which are predicted to 
decay into both J/^ + p and J/^ + uj, as indeed seems to 
be the case [2(| for X(3872). A precise measurement of 
the branching ratios can provide a determination of sin 9 
and therefore a precise prediction of the mass difference. 

We analyze, in this context, the process in which the 
light vector meson from X decay goes into a lepton pair, 
with p — uj interference, which allows to distinguish be- 
tween the two states Xh and Xi. 

Finally, we analyze the non-leptonic decay amplitudes 
B — > KX , for both B + and B° , restricting for simplicity 
to zero mixing. 

From the limit to the width of V(3872) as observed by 
Belle [Hj, we infer that only one particle should domi- 
nate the final state of B + , either X u or X4. The AI = 
rule of the weak transition implies then that B° decay 
is dominated by the other state, Xd or X u : a precise 
measurement of the X mass in B + and B a decay should 
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reveal the mass difference given in @). The observation 
of the decays X — ► J/^P + e + e~, mentioned in the pre- 
vious paragraph, would allow an independent check of 
which particle is which in B + and B° decays. We derive 
also bounds for the production of the charged states 
in B decays, which are close to, but not in conflict with 
the negative results published by BaBar [27j ■ 

It is hardly necessary to remark that our scheme is al- 
ternative to the D — D* molecule picture proposed for the 
AT (3872) [2^|. Albeit in some case one gets similar predic- 
tions (like isospin breaking decays) the particle content 
and the pattern of predictions is quite different, in a way 
that we believe can be put to a test in the near distant 
future. 

The plan of the paper is the following. We discuss in 
Sect. II spin-spin interactions in the constituent model 
and in Sects. Ill and IV the spectrum of the hidden 
and open charm states. Sect. V is devoted to isospin 
breaking, in Sect. VI we discuss the X decays. The 
production of X states in non-leptonic decays of £> +, ° 
is discussed in Sect. VII. We present our conclusions in 
Sect. VIII. 
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490 


1670 


362 


546 


1721 



TABLE I: Constituent quark masses derived from the L = 
mesons (first row) or from the L = baryons (second row). 
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315 


195 


121* 


70 


72 


59 


(Kij)omimj(Ge'Vy 


0.029 


0.029 




0.036 


0.059 


0.16 



TABLE II: Spin-spin couplings for quark-antiquark pairs in 
color singlet from the hyperfine splittings of L — mesons 
(first row). The values in the second row show the approx- 
imate scaling of the couplings with inverse masses (masses 
from meson spectrum) . *The ss coupling which is not exper- 
imentally accessible, is obtained by rescaling the sq one by 
the factor m q /m a . 
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II. CONSTITUENT QUARKS AND SPIN-SPIN 
INTERACTIONS 

In its simplest terms, the constituent quark model |2j, 
l3(i| derives hadron masses from three ingredients: quark 
composition, constituent quark masses and spin-spin in- 
teractions. The Hamiltonian is: 



H = rrii + 2/ty (Si ■ Sj) 



(3) 



and the sum runs over the hadron constituents. The 
coefficients Kij depend on the flavor of the constituents 
2, j and on the particular color state of the pair. 

It is not at all clear how this simple Ansatz can be 
derived from the basic QCD interaction, in particular 
how comes that the effect of the spin-independent color 
forces, responsible for quark confinement, can be sum- 
marized additively in the constituent masses. However, 
it is a fact that Eq. © describes well the spectrum of 
mesons and baryons, with approximately the same val- 
ues of the parameters for different situations. The spin- 
spin interaction coefficients scale more or less as expected 
with constituent masses and, when compared in different 
color states, with the values of the color Casimir coeffi- 
cients derived from one-gluon exchange (as we shall see, 
this is less accurate). Be as it may, we shall accept the 
simple Hamiltonian of Eq. ©. The rest of the Section 
is devoted to determining the parameters from the me- 
son and baryon masses. We summarize the (well known) 
mass formulae for the sq and sq pairs (throughout the 
paper: q — u,d) and give a summary of the parameters 
in Tables I to III. 

Applied to the L — mesons, K and K* , Eq. © 



M 



J(J + 1) 



(4) 



Adding the analogous equations for it — p, D — D*, 
D s — D* , we find four relations for the constituent masses 
and the values of the four couplings, as reported in Ta- 
bles I and II. There is one consistency condition for the 
constituent masses, which can be written as: 



(m c + m q ) D + (m s + m q ) K 
(m c + m s ) Ds = 2076 MeV 



(2m q )ir 



2157 MeV 
(5) 



This relation is representative of the inaccuracy of the 
model. Spin-spin interactions scale as expected, like the 
inverse product of the masses of participating quarks, a 
most remarkable feature. 

Adding the J — r\ c complex, we obtain the cc cou- 
pling, also reported in Table II, and a considerably 
smaller constituent charm mass: 



(mjj/y = 1534 MeV 



(6) 



The parameters from the J/^ system deviate appreciably 
from the rest, a not unexpected feature since the charmo- 
nium wave function is determined by the charmed quark 





qq 


sq 


cq 


cs 


( K y)g (MeV) 


103 


64 


22 


25 


(Kij)3m i mj(GeV)' > 


0.014 


0.013 


0.014 


0.024 



TABLE III: Spin-spin couplings for quark-quark pairs in color 
3 state from L — baryons. One gluon exchange implies 
(nij)a = l/2(Kij)o. The values in the second row, show 
the approximate scaling of the couplings with inverse masses 
(masses from the baryon spectrum). 
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mass and therefore is considerably different from those 
of the mixed flavor mesons, which are determined by the 
light quark masses. 

Baryon masses allow us to obtain quark-quark spin 
interaction in a color antitriplet state. We consider the 
uds states, A ("good" diquark, S = 0), S and Y* ("bad" 
diquark, 5 — 1). One finds: 



M(S, J) = 2m q + m s + (K qq )^ 

+ ( K qs)3 



5(5+1) - 



J(J+l)-S(S+l)-- 



(7) 



We can write similar equations for P — A + , involving 
only (re M )g, and for A c , £ c , £*, involving (/s 35 )g and 
(Kqc)a- We find three determinations of (n qq )s, the values 
of {K qs )3 and (K gc )gj as wen as a new determination of the 
three constituent masses. The new information is given 
in Tables I and III. The consistency conditions read: 



(« w ) g (P-A) 

(«?g)s( A ~ S ~ 
{ K qq)3~{A-c ~ £c 



= 97 MeV 

S*) = 103 MeV 

- K) = 107 MeV 



(8) 



For completeness, we consider also the three H c states, 
involving (K qs )<i again and (K sc )g. From the masses, we 
find: 



(K gs ) s (E c ) = 78 MeV 
(k sc )s(S c ) = 25 MeV 



(9) 



The overall agreement is quite satisfactory, in particu- 
lar for the spin-spin couplings. The decreasing strength 
with increasing mass is evident, until we go to cs(s) or cc 
states which may have considerable distortions in their 
wave functions. 

For our purposes, however, we need to consider fur- 
ther couplings, which refer to the quark-antiquark inter- 
actions to which we have not yet experimental access. 
Inside our states, these pairs are in a superposition of 
color singlet and color octet. Omitting spinor and space 
time variables, we write: 



[cq] [cq] 



^abc , 



e"^e-ab'c'{c b q c ){c b q c ) 
(c b q c )(c b q c )-(c b q c )(c c q b ) (10) 



Color indices in the last term of Eq. (|10fl can be rear- 
ranged with the use of the familiar color Fierz-identities: 



It is not difficult from Eq. I|ll|) to see that the probabil- 
ity to find a particular qq pair in color octet is twice the 
probability of the color singlet, so that (the same holds 
for the other flavors as well): 



1 2 

K c5 {[cq][cq]) = -(Kcc)o + g( K ce)8 (12) 



Of course, we do not know (k cS )s- We resort to the 
rule derived from one-gluon exchange: 

(kcc)x = const. [C^(X) - C^(3) - C< a >(3)] (13) 

where C*- 2 -* (X) is the value of the quadratic Casimir op- 
erator in the representation X: C (2) (X) = 0,3,4/3,4/3 
for X = 0, 8, 3, 3. Eqs. O and {T3J| give, in conclusion: 



K c e = K cS ([cq][cq]) = -(K c g)o 



(14) 



We apply the previous results to determine the con- 
stituent mass of light diquarks, considering explicitly the 
case of the ao(980): 



a (980) = [sg] s=0 [sg]s=o 
We write the Hamiltonian according to: 



(15) 



H = 2m [cs] + 2(n sq U(S s ■ S q ) + (S s ■ S q ,)} 

+ 2K sq [(S s ■ S q >) + (S s ■ S q )} 

+ 2k sS (S s -S s ) (16) 

The state given in Eq. I|15fl is not an eigenstate of this 
Hamiltonian, which is diagonalized only within the states 
with different diquark spin composition, see Sect III be- 
low. However, the latter could as well not exist, so we 
content ourselves with the mean value: 

(a \H\a ) = 984 MeV = 2m [sq] - 3(K sg ) g (17) 

and, using the value in Table III, we find: 



m [sq] = 590 MeV 



(18) 



to put into evidence the state of color of the cc pair: 

[<*] [eg] - l(c b q c )(c b q c ) - l -(c\ A c)(q\ A q) (11) 



For <t(480), with (K qq )^, we find 



m M = 395 MeV 



(19) 



Light diquarks constituent are not at all much heavier 
than constituent quarks. 
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III. THE SPECTRUM OF [cq][cq] STATES 



States can be conveniently classified in terms of the 
diquark and antidiquark spin, S cq , S 5 q', total angular 
momentum, J, parity, P, and charge conjugation, C. We 
have the following states. 

i. Two states with J PC = 0++: 



|0++) = |0 C9 ,0 er ; J = 0); 

|0++') = |i C9 ,W;J = 0) 

ii. Three states with J = 1 and positive parity: 

|A) = |O cg ,V;J=l); 
\B) = |l C3 ,0c r ;J=l); 
|C) = |l c , J l er ;J=l) 



(20) 



(21) 



Under charge conjugation, \A) and \B) interchange 
while |C) is odd. Thus the 1 + complex contains one 
C-even and two C-odd states: 



II 



++\ 



1 



V2 



(\A) + \B)): 



\l+-) = —(\A)-\B)); 



r 



= \C) 



(22) 



One can analyze these states in terms of the states with 
definite values for the spin of cc and qq' . The state with 
both spins equal to zero cannot appear, because J = 1; 
among the others, the only one with C = + is that with 
both spins equal to one. Thus, the state |1 ++ ) in Eq. 122(1 
has a definite value of the cc spin, S C c — 1 ■ 
iii. One state with J PC = 2++: 



)++\ - 



|lcgj 1-cg'i J 2) 



(23) 



The 2 ++ state has also S cS = 1. 

Next, we consider the Hamiltonian, which is the same 
as m Eq. (Unj with s -> c: 

H = 2TO[ C? ] + 2(/C C q)g[(5 c • S q ) + (Sc ■ Sq')] 
(Sq-Sq,) 
+ 2K C q[(S c -Sq>) + (Sc-Sq)] 

+ 2k cS (S c -S 5 ) (24) 

The Hamiltonian is diagonal on the 1 ++ and 2 ++ 
states, with eigenvalues: 



M(l ++ ) = 2m [cq] - (/€ c q)3 + -Kqq 



1 

2 



(25) 



M(2 ++ ) = 2m [cq] + (K cg )3 + \^qq + K C q + ^K c5 (26) 

A tedious but straightforward calculation (see Ap- 
pendix) leads to two, 2x2, matrices for the other states. 



M(0++) = 



2 



-3(k C9 )c 



2 (,^qq ~^ ^cc ^^cq) 



cq)3 



2ft C q) 



M(l+-) 



~i. K cq)z + 



3 T n-CIJ 



2 



(« 



rL qq 



The state 1 ++ is an almost perfect candidate to explain 
the properties of A (3872): 

- it is expected to be narrow, like all diquark- 
antidiquark systems below the baryon-antibaryon thresh- 
old; 

- the unnatural spin-parity forbids the decay in D — D, 
which is not observed; 

- it can decay in the observed channels J/^ + 
light vector meson, with conservation of the spin of the 
heavy flavor system; 

- it decays into both p and u>, due to isospin breaking 
in its wave function (Sect. IV). 

How narrow is narrow we shall consider in Sect. VI. 
For the moment, we identify the 1 ++ with the A(3872) 
and proceed to compute the spectrum via the couplings of 
Table III, directly, and those of Table II, scaled according 
to Eqs. dJ, |I2> and (H3J). By Eq. the diquark 

constituent mass is fixed to be: 



cq] 



1933 MeV 



(27) 



We report in Fig. 1 the full spectrum computed nu- 
merically. The energy levels have an error which is dif- 
ficult to quantify at the moment, maybe in the order of 
10 — 20 MeV. A few observations are in order, 
i. From Eqs. l(23)l and j|2l)|l we read: 



M(2 



M(l 



2[(K cq )- s 



3952 MeV (28) 



This places the 2 ++ close to the recently observed 
resonance at 3940 MeV. Note that the coupling K cq is well 
determined and the other, n cq , could easily be smaller 
than we estimate with the color factor. The identification 
of the 2++ with the A (3940) is quite attractive. The 
2 ++ can decay in J/^+light vector meson respecting the 
conservation of the heavy flavor spin and also in D — D. 
The decay A(3940) J/*+w is seen by Belle, the D-D 
decay should be searched for, but it could be somewhat 
suppressed by the decay in D-wave. 
ii. Of the two ++ states, one is below the D — D thresh- 
old. It can decay in r] c iT or r\ c r\ or multihadron states. 
The other should be seen to decay in D — D. There are 
no candidates, at present, for the ++ states. The same 
holds for the two 1^ states. Allowed decays of the latter 
are J/* + 7r(r?), r) c + p(oj). 
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FIG. 1: The full spectrum of the X particles. 



IV. THE [cq][sq'] STATES 

We extend the calculation of the previous Section to 
the states [cq] [sq'] , leaving aside the issue whether they 
can bind or not. The appropriate Hamiltonian is: 

H = rri[ cq ] + m,[ aq ] + 2(K cq )^(S c ■ S q ) + 

-\-2{K S q)^(Sg ' Sq') + 

-\-2Kqq(Sq • Sq') + 2K C q[S c ■ Sq') + 2n S q(Sg " Sq) + 

+2k c - s {S c ■ S s ) (29) 



The angular momentum composition of the multiplet 
is, of course, the same as the previous one except that 
this set of states is not invariant under C-conjugation 
and the J p = 1 + states form an irreducible complex. 
The energy levels are given by the following formulae. 



M(2 H 



™>[cq] + m [sq] + -[(«c,)a + H^sqh} + 



— K 

2 



2 Kc 



(30) 



M(l+) u 
M(l+) 12 

M(l+)i3 

M(l+) 22 
M(l+) 23 
M(l+) 33 



[ 3(K C q) 3 + {K S q) 3 ]/2 
{^qq ^cq ^sq ^cs)/2 

i^qq K cq K qs K cs)/ V^2 

[( K cq)a — 3( K sq)3,]/2 

{~ K qq ~ K cq + K <jg + K cg)/ V2 
[(^09)3 + { K sq)3 Kqq Kqg — 



(31) 



Kcs]/2 



M(0 + )i 2 



2 [t K '-'/ ' 
2 



M(0+) u = - ; [Ws + Wil 

~( K qq K cg K sg + K cs) (32) 



1 



M(0+) 22 = -[( Kc? )3 + ( Ksg )3]- 

[ftqq H~ ^cg ~f~ H~ ^cs) 

The spectrum is reported in Fig. 2. 
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FIG. 2: The predicted spectrum of the particles with open 
charm. 



The particle claimed by the SELEX Collaboration 
fits quite naturally in it as the 2 + member of the multi- 
plet. Note the change of attribution, with respect to the 
+ assignment suggested previously which makes 
the A (2632) compatible with the diquark constituent 
masses found in Eqs. I|18|) and l|27|l . without changing 
the results presented there. In addition, we associate 
tentatively the lowest + and one of the lowest lying 1 + 
with D s (2317) and D s (2457) [l^, respectively. This is 
compatible with the observed decays: 



£ s (2317) 
D s (2457) 



■ D sl ^- (D s )*7r 



(33) 



A four quark interpretation of the D s particles has 
been advanced in Rcf. 20], while the cs interpretation 
is pursued in [2l|. We find quite suggestive that by as- 
signing the A (3872) to its natural, 1 ++ , level and using 
reasonable values of the spin-spin couplings we are able 
to fit other four particles, which could be at odds with 
the conventional quark-antiquark interpretation. 



V. ISOSPIN BREAKING 

We consider in this Section the finer structure of the 
A (3872). In particular, we consider the neutral states 
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with the composition given in Eq. (JTJ. Physical states 
could be expected to fall in isospin multiplets with I = 
1,0: 

f ce =(X u + X d )/V2; 

a cs = (X u - X d )/V2 (34) 

The two states in Eq. Q are mixed by self-energy di- 
agrams whereby a light quark pair transforms into an- 
other one by annihilation into intermediate gluons. In 
the basis Eq. Q annihilation diagrams contribute equally 
to all entries of the mass matrix. The contribution of 
quark masses, on the other hand, is diagonal in the basis 
Eq. QJ. The resulting 2x2 matrix is: 

/ 2m„ + S S \ 
^ 6 2m d + 8 J 

5 being the contribution from annihilation graphs. The 
matrix with all equal entries 6, admits the states in <|34|) 
as eigenvectors, with split masses. 

At the mass scale determined by the cc pair we ex- 
pect annihilation diagrams to be small, as indicated 
by the very small J/^ width. Thus, mass eigenvectors 
should align to the quark mass basis. For the strange 
quark, this happens already at the mass scale of the vec- 
tor mesons, <f> and u). The other case is provided by the 
scalars a(980) and / (980), which are quite degenerate in 
mass. The upper bound [2<| \AM\ < 10 MeV indicates 
that annihilation contributions are, at best, at the level 
of the normal isospin breaking mass differences, suggest- 
ing a sizeable deviation from the isospin basis. At the 
A(3872) scale, we expect the u — d quark mass difference 
to dominate and the mass eigenstates to coincide with 
the states in (QJ to a rather good extent. 

A numerical estimate of the mass difference is obtained 
as follows. The up and down quark mass difference is 
determined by the pseudoscalar meson spectrum |25j| . 
after separating its contribution from the background 
of second order electromagnetic corrections due to one- 
photon exchange. The so-called Dashen's theorem |3l| 
states that one-photon exchange does not contribute to 
the isospin breaking, U-spin singlet combination of Kaon 
and pion mass differences, which is therefore given by the 
quark mass difference: 

(K+ - K°) - (tt+ - tt ) = C(m u - m d ) = 

= -5.3 • 10~ 3 (GeV) 2 (35) 

where particle symbols stand for squared masses. Com- 
bining with the equations for the pion and Kaon masses 
in terms of quark masses, and assuming m s — 150 MeV 
from the baryon mass differences, one finds: 

(m„ - m d ) = 3.3 MeV (36) 

Before translating Eq. (|36[1 in hadron mass- differences, 
one must control the one-photon exchange contributions. 
We can divide the e.m. corrections in (i) corrections 



on the same constituent quark line, (ii) photon cross- 
ing from one to another quark line. In the constituent 
quark model, the first correction goes into a renormaliza- 
tion of the constituent mass while the second one adds 
to the spin-spin interaction. A control case is that of the 
charmed mesons: 

M(D+) - M(D°) = 4.78 ± 0.1 MeV 
M(D*+) - M(D*°) = 3.3 ± 0.7 MeV (37) 

Note that the result in the first line is larger than the 
mass difference in <|36ll : the spin-spin interaction in total 
spin zero is repulsive (attractive ) for non- vanishing (van- 
ishing) total charge. The correction to the constituent 
masses is obtained by averaging over the two spin multi- 
plets, see Eq. |gjh 

_ 3(£>*+ - D*°) + {D+ - D°) _ 

y^d THu) const ~~ ^ 

= 3.7 ±0.7 MeV (38) 

The agreement with Eq. I|36|l is better than for individual 
mass differences. Unfortunately, at the moment, we do 
not have enough masses to determine and subtract the 
spin-spin e.m. interaction for the [eg] [eg'] multiplet. We 
take the example as suggestive that neglecting photon 
exchange may introduce an error of, perhaps, ~ 30%. 

Non-negligible gluon annihilation diagrams mix X u 
and X d and increase the mass difference. Writing: 

A"iow = cos 9X U + sin 8X d 

Xhigh = — sin 9X U + cos 9X d (39) 

we get the result already stated in the Introduction: 

M{X h ) - M(X t ) = 2(m d - m„)/cos(20) = 
= (7±2)/cos(20) MeV 

The mixing angle can be determined from AM as well 
as from the ratio of the decay rates in J/^+cu and J/*f> + 
p, as we shall see in the next Section. It goes without 
saying that the same considerations can be applied to all 
states in Figs. 1,2. 

In conclusion, we predict close to maximal isospin 
breaking in the wave function and correspondingly in the 
hadronic decays of V(3872). 

Isospin violation in the wave function is also predicted 
by the DD* molecule scheme |28|, with the AT (3872) 
being essentially (D°D*°) + (D°D*°). However, in our 
scheme we have two states rather than one, separated by 
the mass difference {J2J, and quite a richer phenomenol- 
ogy. 

VI. THE V(3872) DECAY WIDTH 

A pair of color-singlet mesons cannot be obtained by 
cutting the strings that join quarks and antiquarks in the 
H shaped [qq] [qq] states. The baryonium picture suggests 



8 



that the two-meson decays of the latter go via intermedi- 
ate baryon-antibaryon states of high mass. This implies 
basically narrow widths. 

The X(3872) is expected to be particularly narrow for 
several additional reasons. 

(i) Unnatural spin-parity forbids decays into DD; 

(ii) the channel DD* is below threshold; 

(iii) decay in rj c + mesons is forbidden by heavy flavor 
spin conservations |l6j . 

Of the charmonium channels, the only available ones 
are J/^f + 2ir and J/*f> + 3ir, dominated by p° and oj, 
respectively. Each mass eigenstate decays simultaneously 
in the two channels, due to isospin breaking in the wave 
function. 

We describe the decay by a single switch amplitude, 
associated to the process: 



CUsCUl 



(cc)o(uu) 



(40) 



where subscripts indicate color configurations. 

We further write the invariant three-meson coupling 
for X u according to: 



= g v M x (XAiP)-V 



(41) 



where P p and M x are the decaying particle momentum 
and mass. To estimate the value of gy, we compare with 
the similar couplings for the light scalar mesons, deter- 
mined by one dimensionful constant, A ~ 2.6 GeV. An 
admittedly bold guess, to obtain the order of magnitude, 
is: 



g v M x 



V2 



(42) 



By dominating the 2tt(3tt) decay with p(ui) exchange 
in the narrow-width approximation, we find: 



dT(X l ^^j + f) 2x lyV \A\ 2 B (v . 



•I) 



ds 



8irM x 



M V T V 



p(s) 



(s - M 2 )2 + (MvTvY 



(43) 



with / = 7r + 7r (tt + tt 7t°) for V = p(co), s the invariant 
mass-squared of the pions, and p the decay momentum: 



p(s) 



y/X{M x ,M^M v ) ^ 
2M X 5 

\4 i C\f ^4 , fur \4 n /n,r n,r ^2 



A = (M x ) 4 + (M^) + (My) - 2(M X M 4 ,) 2 
-2(M x M v f - 2(M il M v ) 2 (44) 



The coefficient xi y is: 



XLV 



(cos9±smBf 



(45) 



for V = oj(p). Similar equations hold for the higher mass 
state, Xh, with the appropriate substitutions. 



By numerical integration, we then find: 



(p) P = 



* / J(2m„r ( s ~ M 2 ) 2 + (MpTp) 2 
= 126 MeV; 

(p) u = 22 MeV (46) 



and: 



= 2x hp ■ 2.3 MeV; 

T(X l ^J/^ + ^-^) = ^^(p) u 



SwM 



x 



2xi „ ■ 0.4 MeV 



(47) 



We anticipate small widths, comparable to the resolu- 
tion of Belle and BaBar. However, given the mass differ- 
ence in l|2"jl. one would expect to observe either two peaks 
or one unresolved structure, broader than the stated ex- 
perimental resolution, 4.5 MeV. Taking Belle data at face 
value, we conclude that only one of the two neutral states 
is produced appreciably in B + decay (this will be dis- 
cussed in the next Section). Assuming this to be the 
case, we can get some information on the mixing angle 
from the observed ratio of Zir to 2ir decay rates. We get 
from Eq. |g7J|: 



n^)y Xl 
n^)) Xh 



(cosfl + sinfl) 2 (pj) 
(cos 9 — sin 9) 2 (p p ) 



(cos 9 



(cos 9 + sin 9) 2 (p p ) 



(48) 



750 MeV to u> decay and divides by the total number of 
observed 2ir events. They find: 



/r(37r; 



= 0.8 ± 0.3 ata t ± 0.1 



syst 



(49) 



Vr(27r) /Bcllc 

The central value is compatible with Eq. 148|) for: 

9 ~ ±20° (50) 

for Xi or Xh, respectively. Assuming that there are no 
other significant decay modes, the corresponding widths 
and branching fractions for the particle seen in B + decay 
are: 



T = 1.6 MeV (3.7 MeV) 
B(2tt) = 0.61 (0.95) 



(51) 



where we have listed in parenthesis the properties of the 
particle not seen in B + decay. The mass difference of the 
two states is: 



M(X h ) - M(Xi) = (8 ± 3) MeV 



(52) 
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We give also the corresponding predictions for the 
charged state X + , which decays via p-exchange only: 



T(X+ -> J/V> + = ^r(p) P = 4.6 MeV; 



8nM 



x 



B(X + -> J/7/. + 7T + 7r°) ~ 1 



(53) 



The value of the mixing angle in (|50|l is perhaps on the 
high side but still compatible with the general picture. 
More precise data are clearly needed. 

We close this Section by considering the leptonic de- 
cays: 



X(3872) -> J/* 



(54) 



The lepton pair originates from the coherent superpo- 
sition of p and u> produced in the decay of the X . Thus, 
the branching ratio can distinguish between Xi and Xh, 
supplementing the measurement of the mass. For sim- 
plicity we give the result for the case of vanishing mixing. 
A simple calculation gives: 



dT(X -> ip + 



e+e~) 



ds 

\A\ 2 B(p^e+e-)M p T p 
87rM| 7T 



± 



p(s) 



1/3 



(s - A/2) + i{M p T p ) (s - A/2) + »(M w r w ) 

We have assumed the quark-model ratio for the leptonic 
amplitudes of p and u> and used the narrow width approx- 
imation. The sign ± applies to X u and X d , respectively. 
Combining with Eq. (|47|l . with 6 = 0, we find: 



B(X U 
B(X d 



0.8- 1CT 4 
0.3- 1CT 4 



(55) 



VII. PRODUCTION OF [cq][cq] STATES IN B 
NON-LEPTONIC DECAYS 

B + and B° decays produce superpositions of the two 
neutral states, Eq. Q as well as the charged states: 



X+ = [cu][dc]; X- = [cd\[uc] 



(56) 



For simplicity, we shall restrict to vanishing mixing. 

We consider first the B + decay amplitudes for the al- 
lowed decay: 



(bu) ^c + c+ s + u+(u + uord + d) (57) 



One additional pair is included in the final state, created 
from the vacuum by the strong interaction. 

If we want a K in the final state, the s must combine 
either with the spectator quark, u, to give a K + (am- 
plitude A\) or with one quark from the additional pair, 
to give either K + or K°, (amplitude A2). Thus we have 
two independent amplitudes: 



B+: 



A(K+X U ) 
A{K+X d ) 

A(K S X+) 



A^ 

An 
M 

' V2 



A 2 



(58) 



For B decays we have simply to exchange u with d 
and K+ with K°, to get: 
B°: 



A(K s X d ) 



Ax + A 2 



V2 



A(K S X U ) = A ; 
A(K+X-) = A 2 



(59) 



We note, in passing, that these relations follow also from 
the AJ = rule obeyed by the weak transition. 

We noted already that the mass difference given in 
Eq. (J2J is larger than the apparent width of the X(3872) 
peak seen by Belle Thus, only one of the two neu- 
tral states is produced appreciably in B + decay. Orien- 
tatively, we shall assume that: 



T(K+X 



(u or d) 



>AT(K+X idolu) ) 



(60) 



Eq. H60JI implies some bound to the production of the 
charged states, X , not observed thus far. With three 
amplitudes and two parameters, Eq. (|58|l gives rise to the 
triangle inequality: 

\A(K+X U )\ + \A(K+X d )\ > V2\A(K s X+)\ 
>\\A(K+X u )\-\A(K+X d )\\ (61) 

We are interested in the lower bound to the rate of X + , 
which, due to Eq. (JBUJl or Eq. (JSTJ, is: 



1 



\A(K s X+)\ > —\A(K+X q )\ 



(62) 



with q — u or d, according to which is the dominant 
decay product. 

Eq. H6Uf) has two solutions: 
B + — > K + X„ dominant: 



Ax ~ A 2 



or: 
B+ 



K + X d dominant: 



Ax 



(63) 



(64) 



We consider now B° decays, Eq. I|59|l . It is immediate 
to see that if X u dominates B + decays, X d dominates 
B° decays and viceversa. Thus we are led to predict that 
the X particle in B + and B decays are different, with 
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a mass difference given by Eq. or Q52JI. In addition, 
from the corresponding triangle inequality, we find: 



T(K+X-) > -T{KsX q ) 



(65) 



with q = u or d, whichever particle dominates B° decays. 
Relations i|62l65|) remain unchanged if one considers an 
equal mixture of B° and B° decays and adds K + and 
K~~ events. 

To conclude, we give explicitly the lower bounds to the 
production of X ± in B + and B° decays: 



B(B+ -» K S X+) ■ B(X+ -» J/fr + tt+tt 
6(5+ - K+X l/h ) ■ B{X l/h -» J/* + 7T+7T 



> 0.2 



> 0.53 



B(B° -> K s X h/l ) ■ B(X h/l - J/* + tt+tt-) 
to be compared with the upper limit given by BaBar |27| : 



R + < 0.8 



with large errors. 



VIII. CONCLUSIONS 

The diquark-antidiquark structure explains well the 
properties of the X(3872): spin-parity, narrow width, 
simultaneous decay into channels with different isospin. 
Taking the X(3872) as input, we have derived a spectrum 
which is able to explain spin-parity and decay properties 
of few other particles that could be at odds with a qq pic- 
ture: X(3940), the previously discovered Z? s (2317) and 
D s (2457) and the X(2632) claimed by SELEX. 

Isospin breaking in the wave function and in strong 
decays of these states is a distinctive consequence of 
the asymptotic freedom of QCD, much in the same way 
as narrow widths for heavy quarkonia. Also, all these 
states have to be doublets, unlike the case of the D — D* 
molecules, with typical mass splittings given by twice the 
down-up quark mass difference. The two different states 
of X(3872) should appear in B + and B° decays, respec- 
tively. 

The crucial test of the scheme, of course, will be the 
observation of the charged or doubly charged partners of 
the X particles and, more generally, the observation of 
heavy states with really exotic quantum numbers. We 
have derived rather strict bounds for the production of 
X + (3872) in B decays, close to the present limits so that 
a meaningful test may be expected in the near future. 
The existence of exotic states at low-energy is also a 
pressing issue. 

The indications derived from the properties of 
a//(980) and X(3872) seem to us very compelling, so 



as to warrant a thorough experimental and theoretical 
investigation. 
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APPENDIX 

We give here a simple derivation of the matrix elements 
of the spin-spin Hamiltonian over the qq — qq states. We 
consider states which can be written schematically as: 

\S lcu] ,S [sa] ; J) = |T, T'; J) = {c a T ab u b )(c c T' cd u d ) (66) 

Sjcu] and S[cu] are the total diquark and antidiquark 
spin and J the total angular momentum. Individual spins 
are represented by 2 x 2 matrices, P, with: 



r = 71' r = V = 2 a2a 



(67) 



for spin and 1, respectively. The matrices T are nor- 
malized so that: 

Tr[(r Q )tr /3 ] = S af3 (68) 

Spin operators are defined according to: 



Since: 



s, l |r) = |ria) ; s c |r) = |ia T r) 



a a 2 = —cr 2 cr 



(69) 



(70) 



we recover the expected formulae for the total spin oper- 
ator: 

(s„ + s c )|r°) = o 

[(S u ) i + (S c ) i }\P)=ie^ k \T k ) (71) 

We find, also: 



(0|S U |1) = -(0|S C |1) = -; 

(l|S„|l) = +(l|S c |l) = i(l|(S„ + S c )|l) (72) 

We can now compute the matrix elements of products 
of spin operators. We have two cases. 
1. Same diquark, e.g. S u ■ S c . This operator is just a 
combination of Casimir operators: 

2(S U • S c ) = (S cu ) 2 - (S c ) 2 - (S u ) 2 (73) 
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and is diagonal in the basis. 

2. Different diquarks, e.g. S u ■ S a - We consider as an 
example the J = states, represented by (summation of 
repeated indices understood): 

1 0, 0; 0) = l - (a 2 ) ® (a 2 ) ; 1 1 , 1 ; 0) = ^= (a 2 a 1 ) ® (<r 2 <r 4 ) 

(74) 

Using the basic definitions, we have: 

2(S U • S u )|0, 0; 0) = Ifaj) ® (a 2 a l ) = ^|1, 1; 0); 

2(S U • S u )|l, 1; 0) = -L(<7 2 aV) ® (ffjffV) = 
4V3 



v/3 



0,0; 0> - |1,1;0) 



(75) 



In conclusion, on the states |0, 0; 0) and |1, 1; 0) we obtain 
the matrices: 



2(S U • S c ) — 



-3/2 
1/2 



2(S M • S u ) 



y/3/2 
y/3/2 -1 



(76) 



Using the relations, Eqs. (|72|l . we derive from l|7ti|) the 
representatives of the other spin-spin operators, to ob- 
tain, e.g. the mass matrix given in Sect. III. 

We conclude by giving the tensor basis for the J = 1 
states. 



\A) = 


|0,i;i; 


\B) = 


|1,0;1 


\C) = 


|i,i;i 



2V2 
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